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Abstract 

O ; 

A system of charged bosons at finite temperature and chemical potential is studied 
| in a general-relativistic framework. We assume that the boson fields interact only 

■ gravitationally. At sufficiently low temperature the system exists in two phases: the gas 

and the condensate. By studying the condensation process numerically we determine 
the critical temperature T c at which the condensate emerges. As the temperature 
decreases, the system eventually settles down in the ground state of a cold boson star. 
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1 Introduction 



Bose-Einstein condensation is an extensively studied subject of interest to almost all branch- 
es of physics. It is probably the simplest example of critical phenomena since it exists 
theoretically even in noninteracting systems [|l], ||. In an interacting theory containing 
complex scalar fields, Bose-Einstein condensation becomes a highly nontrivial phenomenon. 
In self-interacting scalar theories there are cases, such as 4 theory, where a homogeneous 
condensate is a stable ground state. In some cases, the Bose-Einstein condensate at zero 
temperature exists in the form of stable nontopological solitons known as Q-balls [|], |J. An- 
other interesting example is a scalar field coupled to gravity. The ground state of a condensed 
cloud of charged bosons of mass m, interacting only gravitationally and having a total mass 
M below a certain limit of the order MpJ m, is a stable spherically symmetric configuration 
H which is usually referred to as mini-soliton star [0] or boson star ||. The gravitational 
collapse of boson stars is prevented by Heisenberg's uncertainty principle. Seidel and Suen 
have demonstrated that boson stars may be formed through a dissipationless mechanism, 
called gravitational cooling ||. Boson stars have recently attracted much attention as they 
may well be candidates for nonbaryonic dark matter [|HJ . 

In this paper we study the thermodynamics of boson stars and its relation to the Bose- 
Einstein condensation in the framework of general relativity. Some aspects of boson stars 



at finite temperature have been studied in Newtonian gravity ||11||. Our approach is similar 
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in spirit to the analysis of Q-ball thermodynamics by Laine and Shaposhnikov fl2" |. We will 
consider a canonical ensemble of charged bosons and study its properties near the critical 
temperature at which the Bose-Einstein condensate forms. 

We organize the paper as follows: In Section we discuss the grand canonical ensemble 
of gravitating bosons within a general-relativistic framework. In Section |3| the free energy is 
derived for a canonical system of gravitating bosons. In Section |4] we discuss the equations 
that describe a boson star at finite temperature and their numerical integration. Numerical 
results are presented in Section [5]. Conclusions are drawn in Section |6]. 

2 Grand canonical ensemble 

In this section we derive the grand canonical partition function for a system of self-gravitating 
charged bosons contained in a sphere of large radius R in equilibrium at finite temperature 
T = 1/(3. In a grand canonical ensemble we introduce the chemical potential fi associated 
to the conserved particle number N. The partition function is given by 

Z = Tre-^ H -^ = J [dg}ld$][d$*]e- s z- SKG i (1) 

with the Euclidean actions S g and Skg for the gravitational and the Klein-Gordon fields, 



respectively. The gravitational part may be put in the form |13| , |14 



S g = -T^- f d 4 x^TZ - ±- / d 3 xVh(K - Ko), (2) 
167T Jy Sir JdY 

where 1Z is the curvature scalar, h is the determinant of the induced metric on the boundary, 
and K — Kq is the difference in the trace of the second fundamental form of the boundary 
dY in the metric g and the flat metric. The boundary is a timelike tube of radius R which 
is periodically identified in the imaginary time direction with period (3. Thus, the functional 
integration assumes the periodicity in imaginary time and the asymptotic flatness of the 
metric fields. The matter action is given by 

S KG = J y d 4 x^g(g^d^*d u <S> + m 2 |$| 2 ), (3) 



with the replacement of the Euclidean time derivative |T5| 



where the + or — sign is taken when the derivative acts on <3>* or $, respectively. The path 
integral is taken over asymptotically vanishing fields which are periodic in imaginary time 
r with period (3. The dominant contribution to the path integral comes from metrics and 
matter fields, which are near the classical fields g^, $*, and $. The classical fields extremize 
the action, i.e., they are solutions to the classical field equations. The background metric 
generated by the mass distribution is static, spherically symmetric, and asymptotically flat, 
i.e., 

ds 2 = e^ r) dt 2 - e x(r) dr 2 - r 2 {d$ 2 + sin 2 ?? dip 2 ). (5) 
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Substitution t = —ir converts this into a positive definite metric with the Euclidean signa- 
ture, 

dsl = e u{r) dr 2 + e x{r) dr 2 + r 2 (d$ 2 + sin 2 $ dip 2 ). (6) 
In order to extract the classical contribution, we decompose <£> as 

$0) = <f>(x) +ip{x), (7) 

where is a solution to the Klein-Gordon equation which we will call condensate and ip 
describes quantum and thermal fluctuations around 0. If we neglect quantum fluctuations 
of the metric, the partition function factorizes as 



Z = Z g Z cd |[#][#*K 



■Skg[W 



where Z g = e Sg and Z c & = e Scd represent the saddle-point gravitational and condensate 
contributions, respectively. The condensate action is given by 



S cd = l^x^g^d^d^ + mM 2 ). (9) 



Variation of S c d with respect to <fi* yields a Klein-Gordon equation in which the time deriva- 
tives are replaced by (f|). With the metric (|6|) and the usual ansatz for static solutions 

0(t, r) = ^e--'V(r) = ^e~"\(r) , (10) 

<f>*(t, r) = ^e^'V(r) = , (H) 

we obtain the static Klein-Gordon equation 

d 2 <P _|_ 1 (dv dX ^ 4\ dip ^ A 2 _ v 2 ^ 
dr 2 2 \dr dr r J dr 

where u = to' + \i. Regular, asymptotically vanishing solutions to equation (|T2]) coupled 
with Einstein field equations are solitons which describe the well-known boson stars |7], [8| 
at zero temperature. Our aim is to extend and analyze the corresponding solutions at finite 
temperature. 

Ignoring for the moment the periodicity condition, the ansatz (050) gives the condensate 
contribution to the partition function 

In Z cd = ±/3 jf Sx^g [gVp 2 - g^ip) 2 - m 2 p 2 ] , (13) 

where E is a spacelike hypersurface that contains the condensate. Using this expression we 
find the net number of particles in the condensate 



ldlnZ cd 



j[ d 3 x^g 00 u;p 2 . (14) 
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Alternatively, we define the particle number in a covariant way as 

Ned = J n^dT,^ = J d 3 x y /g^n cd , (15) 

where n c d is the particle-number density in the condensate and g^) = det(^), i, j = 1,2,3. 
Here we have used the fluid four-velocity w M the components of which in the comoving frame 
are 

u, = ^. (16) 



'9oo vS'oo 
Therefore, we identify the particle-number density due to the condensate as 



n cd = vV ^ 2 - (17) 

It may be easily verified that the four-vector j M = n^u^ coincides with the current defined 

as 

j^g^&d^-pd^). (18) 

In order to satisfy the periodicity condition, i.e., <j>{0) = 0(0), we must set uj' = 0. This 
in turn implies uj = fi. In other words, as it was pointed out by Laine and Shaposhnikov 
12]| , if there exist soliton solutions with some values of uj, these solutions are saddle points 



of the Euclidean path integral at fi = uj. In the absence of gravity, a saddle-point solution is 
a homogeneous Bose-Einstein condensate of a free relativistic boson gas. 

Next, we calculate the thermal contribution to Z starting from the expression 

Z th = j[d^][dr]e- s ™ [ ^\ (19) 

with 

S KG = J Y d 4 x^(g^d,rdu^ + m 2 \M 2 ). (20) 

The metric in equilibrium is static, i.e., g^ is independent of r and got = 0. Hence, the 
determinant of g^ u factorizes as g = goog(3)- By making use of the substitution f = r^/g^ 
we obtain _ 

S K g = J d 3 Xy/g\s) df£(f,x), (21) 

where 

C(f,x) = (d f + fiW(d f - + g^d^djip + m 2 \ip\ 2 . (22) 
Here we have defined the local chemical potential and the inverse local temperature as 

P- = /V v^oo! P = PV900- (23) 
These expressions are nothing but the well-known Tolman conditions for chemical and ther- 



mal equilibrium in a gravitational field |16], [17], |18| . Thus, the parameters (5 and \i are 



interpreted as, respectively the values of the inverse temperature and chemical potential 
at infinity fl9fl . In a close neighborhood S of a space point X we introduce local spatial 
coordinates y = fx{ x ) such that 

dxi dxj dyi dyj ' 
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where ij)x{y) = 4 , {fx 1 {u))- The Lagrangian (|22| ) is now locally represented by the Euclidean 
flat-space Lagrangian 

C(r, X,y)=^—+^ x j^—-^ x j+ w — + m \^x I • (25) 
By integrating C(f,X,y) over y we define 

£(r,X) = -^ / d 3 yjr(T,X,y), (26) 

Vo -/So 

where Vo is the volume of S . Using the adiabatic approximation C(f,x) « £(f, X), which 
holds for sufficiently small Vo, we obtain 

S KG = / d 3 Xy/g m (X) ^ df±r d 3 yC(f, X, y). (27) 

Now, the partition function (|T9|) may be factorized as 

Z th = ~[[Z x , (28) 



where 

Z x = J [#][#*]exp[- J d%d%'r(y)K x (y,y')^(y% (29) 

with 



K x (y,y') = V^y/gwW [(d f - fi) 2 + V 2 y )]5^(y - y'). (30) 

Assuming that Vo, although small, is still macroscopic and that the thermodynamic limit 
may be applied to each Z x , the standard functional integration technique |1| gives 



lnZ th = J c^xTrlnfC 1 



= -Jd 3 x V g^ ) J ^[ln(l - e~^) + ln(l - e"^)], (31) 

where E = y/q' 2 + m 2 . This expression may be regarded as a proper volume integral of the 
local partition function lnz(x) 

]nz = - f Apnfl - e"^-^) + ln(l - e -^+A))l, (32) 
J (27r) ,;i 

from which the pressure, energy density, particle-number density, and entropy density may 
be derived in the usual way: 

1 

Pth = p In z, (33) 



9 l 

p th =- ¥ ln, 



(34) 

PP 
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o- = j3 (pth + Pth - ftn th ) . (36) 

These expressions together with (^) yield the well-known parametric representation of the 
equation of state of a relativistic Bose gas in curved space. This equation of state is a special 
case of a more general expression derived from the relativistic kinetic theory fl8fl . 

The gravitational part of the partition function may be calculated from (^) with help of 
Einstein field equations. Using the result of Gibbons and Hawking for the surface term [|13j . 
we obtain 

In Z g = -f3M + J y rf 4 Xv^T °, (37) 

where M is the total mass and is the energy-momentum tensor of the Klein-Gordon field 
averaged with respect to the partition function (|8|). The energy-momentum tensor of the 
Klein-Gordon field can be obtained by variation of the the Euclidean Klein-Gordon action 
with respect to the metric 

= - JL*^G. (38) 

From (0) with we find 

T M „($) = -dfi*d v ® - d v &dfi + g^ig^dp&dv® + m 2 |$| 2 ). (39) 

The averaged energy momentum tensor may be split up into two parts: 

T/ja, =T /JU/ (<p) + T th ^ , (40) 

where the first term on the right-hand side is the classical part which comes from the con- 
densate and the second term represents the thermal and quantum fluctuations. By making 
use of (|10|) and (|ll|) we obtain the condensate contribution to Tq : 

p cd = T °(^) = ifcAV + g U (d i¥ >f + mY). (41) 
It may be shown that the thermal part is of the form which characterizes a perfect fluid: 

?th = (Pth + Pth)UfiU u - PthQuv , (42) 

where the thermal pressure and energy density are given by ([33]) and (0), respectively. 
Putting the condensate (|T3|), the thermal fl31|) , and the gravitational ([37]) contributions 
together, we find the total thermodynamical potential as 

(J) = ~~ In Z = M - jf d 3 x^ [g 00 u;V + p th + p th ] . (43) 
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3 Canonical ensemble 



Now consider a system of self-gravitating charged bosons with the net number of particles 
(number of particles minus number of antiparticles) N contained in a two-dimensional sphere 
of large radius R in equilibrium at finite temperature T — 1//3. In a canonical ensemble, 
instead of the chemical potential /x we fix the particle number which is the sum of the con- 
densate and thermal contributions. Thus, a canonical ensemble is subject to the constraint 

/ d 3 x y /g^(n cd + nth) = N, (44) 

where n cd is given by (|17D and n th by fl35|) with (|32|). 

The free energy of a canonical ensemble may be derived from the grand-canonical parti- 
tion function with the help of the Legendre transform 

F(l3,N)=n(l3,n)+nN. (45) 

The quantity \i in this expression is an implicit function of N and T, such that for given N 
and T the constraint © is satisfied J20|. From © and © with (0) it follows 

F = M-{u-fi)J^ d 3 x^g 00 uj^ 2 - £ d 3 x^7j{p th + p th - /2n th ) . (46) 

If u = /i, the second term on the right-hand side vanishes and the free energy may be 
expressed in the familiar form 

F = M-TS, (47) 
where the total entropy S is defined as a proper volume integral 

S = J au^d^^ = J d 3 x^/g^a, (48) 
over the entropy density a given by (|36|). 



4 Numerical integration 

Given the temperature at infinity T, the radius of the sphere R, and the particle number N, 
we have to solve a set of self-consistency equations consisting of the Klein-Gordon equation 
(|12|) and Einstein field equations: 

87r(p cd +p th )re A + -(e A -l), (49) 
r 

87r(p cd + p th )re A --(e A -l), (50) 
r 

where the thermal pressure and energy density are given by ( j33l) and (|34]) , respectively, the 
condensate density p cd is given by ( |4T| ) and the radial pressure of the condensate by 

p cd = -T r %) = \{gV^ + g%d lV f - mV). (51) 



dv 

dr 
dX 

dr 
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In addition, we impose the constraint which fixes the chemical potential \x at infinity. 
For numerical convenience, let us introduce the substitutions 



,A(r) 
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m 2 x(r) + 1 ' 
1 

1 - 2M(r)/r' 



and new dimensionless parameters 



a 



P 
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to 



P 



n 



m 



(52) 
(53) 

(54) 



where Mpi = yhc/G denotes the Planck mass. Furthermore, we choose appropriate length 
and mass scales such that the length is measured in units of 1/m, and the mass in units of 
Mpj/m. In this way, all physical quantities appearing in our equations become dimensionless 
and equations (0), (f49"D, and (|50| ) now read 



d 2 p ^ ( 2M - 47rr 3 (p - p) \ 1 dp 



dr 



r-2M 



r dr r — 2M. 



XV = o, 



dM A 2 _ 
— — = Arcr p, 
dr 



dx 
dr 



^ = -2(x+l) 



./Vf + 47rr 3 p 
r(r - 2.M) ' 



To these three equations we add 

dAf 
dr 



4vrr 2 (l - 2M/r)~ 1/2 n, 



imposing the particle-number constraint as a condition at the boundary: 

Af(R) = N. 



In equations 



)-(|5q) we have used the abbreviations 
P 



P 



m 2 M 2 x 



(X + 2)^ + 



o , r — 2M ( dip 



dr 



VPth 



V 



V 



m 2 M 2 l 



, r-2M (dp 



r 



+ VPth 



II 



n 



X+ l<p + 77 77, th 



where the dimensionless quantities p t h, Pth, an d 77th are derived from (|32D -(p5] 



~ _ Pth 
Pth — — 7 
m 4 



1 /"°° / 

2^/0 dyy 2 \ l + y 2 ( n + + n - 



(55) 

(56) 
(57) 

(58) 
(59) 

(60) 

(61) 
(62) 

(63) 
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Pth 



Pth 

m 4 



1 



67T 2 JO 



dy 



2tt 2 



dyy 2 



(n + + n. 



71+ 



n. 



n± 



exp 



a( 7v /(l + 2/ 2 )/(x + l)Tl 



1 



(64) 

(65) 
(66) 



In equations (|5l|)-(^6|) the radial coordinate r, the enclosed mass A4, and the particle number 
iV are measured in units of 1/m, m/r], 1/n, respectively, and the scalar field (p in units of 
M P1 . 

It is evident from (|60l) -(|66l) that the thermal contribution to the local quantities is sup- 
pressed by the factor n, of the order of 10~ 38 for a boson mass of the order of 1 GeV. However, 
the contributions of the condensate and the thermal part to the global quantities, such as 
M, N, or F, may be of the same order of magnitude provided the volume of the system is 
sufficiently large. Indeed, the radius R is much larger than the radial extent of the conden- 
sate, and, as we shall shortly demonstrate, the physically reasonable choice of R makes the 
condensate and the thermal contributions comparable in magnitude. We shall also see that if 
we take the limit of vanishing G keeping the central density fixed, we recover the free boson 
gas in which the densities of the condensate and the gas are homogeneous and comparable 
in magnitude for the temperatures above zero and below the critical temperature. 

Equations (|55|)-(|58D should be integrated with boundary conditions dictated by physics. 
To avoid singularities at the origin, we take 



M(0) = 0; W(0) = 



dip 
dr 



0. 



(67) 



The initial value of the scalar field <po = <^(0) at r = is arbitrary but it can be taken to be 
positive without loss of generality. The metric field at the origin xo = x(0) is restricted to 



K xo < Y 



1 



(68) 



because of the requirements that the metric should be positive and that the Bose-Einstein 
distribution (^) should be nonnegative everywhere. However, ip and xo must be simulta- 
neously tuned in order to satisfy the constraint (|59D and to provide the correct asymptotic 
behavior <p(r) —>■ as r — >• oo. The condition that the field (p should vanish asymptotically is 
necessary in order that the condensate be either a soliton or absent. A homogeneous conden- 
sate is excluded for G ^ owing to our assumption that the metric should be asymptotically 
flat. 

Once we have found a solution on the interval (0, R) for arbitrarily fixed parameters a 
and 7 with the correctly chosen initial conditions at r = 0, we can determine the unknown 
physical parameters u, /i, and T from the boundary conditions at r = R. We cut off the 
matter from R to infinity and join the interior solution onto the empty space Schwarzschild 
exterior solution 

4 (i - 

m A \ r J 



x(r) 



- 1. 



(69) 
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Combining this equation with the numerically found interior solution at r = R, we obtain 
u> which then, together with fl54|) , fixes the chemical potential p and the temperature T at 
infinity. 

Consider first the condensate at zero temperature. We know that a zero-node soliton 
solution exist for any particle number N below iV max = 0.653 This maximal value, 
known as the Kaup limit ||, corresponds to the initial values <po = 0.0764, xo — 0.545, and 
the frequency uj = 0.853m. Lower <p lead to lower N and larger uj. In the limit <po — > 0, one 
approaches the Newtonian regime in which uj approaches the limiting value m. It may be 
easily shown that in the Newtonian limit there exist a mass-radius relationship of the form 

M 2 

MR = const— ^, (70) 
m A 

where the radius Rq of the boson star may be conveniently defined, e.g., as the radius of a 
ball of mass M with constant density equal to the central density p(0) = (1 + xo/2)iplm 2 Mp V 
Other definitions are possible [[/J and lead to the same qualitative conclusions. Relation ( [70D 
enables us to analyze the fate of the boson star in the limit G — > 0. In that limit, we keep a 
certain physical quantity fixed, e.g., the particle number N, or the central density p(0). In 
both cases, the radius R$ increases with decreasing G and becomes infinite for G — 0. For 
example, if p(0) is kept fixed, Ro grows as G~ l / A as a consequence of (f70|) . This situation is 
illustrated in Fig. [1] where we plot the number-density profile of the boson star in the ground 
state for several initial values <po- The density profile is in the form of a plateau the length of 
which roughly measures the radius of the boson star. With decreasing <po> which corresponds 
to decreasing G, the plateau widens and becomes more pronounced. In the limit G — > 0, the 
boson star becomes an ordinary homogeneous condensate, with uo = m. 

We now turn to the study of nonzero temperature. As discussed in section ^ thermo- 
dynamic consistency strictly requires 7 = 1 which implies Xo < owing to (|68|) . It may be 
shown that if Xo < and <po > 0, then the function <p(r) is monotonically increasing with 
r. To see this, note that near the origin the first derivative of <p behaves as <p' ~ — <poXo?"/3, 
and hence, <p increases for small r. Furthermore, at every point r > where <p'( r ) = 0' the 
second derivative <p"(r) is negative as a consequence of equation (|55|). This implies that the 
function <p never reaches a maximum. Thus, equation ( |55] ) with (K5|)-(|57|) possesses no soliton 



solution if xo < 0, and at any temperature T > 0, the only physically acceptable solution, 
would be the trivial one <p = 0. It seems that, in contrast to the self-gravitating fermion 



gas [PHI , a self-gravitating boson gas does not possess a nontrivial limiting configuration as 
T — > 0. Moreover, if we hold N fixed and decrease the temperature sticking to 7 = 1, i.e., 
<p = 0, p will increase until we reach a limiting temperature T c for which p(0) = m. Further 
decrease of T with fixed N is no longer possible unless a condensation process takes place. 
If there were no gravity, we would decrease the temperature below T c keeping p = m. As 
a result, the thermal density of particles would decrease and the density of particles in the 
ground state would adequately increase forming a homogeneous Bose-Einstein condensate 
<p = const. In that process, the total number of particles would remain fixed. 

In the presence of gravity, the formation of a homogeneous condensate is not possible 



since <p = const is not a solution to fl55|). However, a condensation process will take place if 
below T c the formation of a soliton is made possible by allowing xo > 0. This implies that 
the chemical potential at infinity p is less than u, although the local chemical potential at 
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the origin p,(0) stays equal to m. The reason for this apparent thermodynamic inconsistency 
lies in the fact that general relativity does not allow configurations with asymptotically 
constant matter density. Of course, in the limit G — > 0, the soliton solution would become 
an ordinary homogeneous condensate, with fi = u = m, as in the T = case discussed above. 
In numerical calculations we shall stick to the condition /2(0) = m below T c . This condition, 
which prevents the Bose-Einstein distribution from becoming negative, fixes 7 = \/xjo + 1 
for Xo > 0, and 7 = 1 for Xo < 0. 

To proceed with the numerics, we have to specify the radius R. Although R is an 
arbitrary parameter, the following physical considerations will give us its preferable order of 
magnitude. First, the size of the system is clearly much larger than the size that it would 
assume at zero temperature, i.e., R must be larger than the radius of a boson star Rq ~ 1/m. 
Second, it should be smaller or at most of the order of the size it would take as a gas at 
high temperatures. To estimate the natural size of the purely gaseous phase, note that if 
the condensate vanishes, the quantities M. and r may be rescaled so that the dimensionless 
field equations remain in the same form as in (0)-(0) with p and p replaced by p th and 
p th , respectively, and with r and M. measured in units of M P i/m 2 and M^/m 2 , respectively. 
Thus, the natural size of a purely gaseous phase is of the order Mpi/m 2 and therefore we 
expect the radius of our system to be in the range 1/m <C R < Mpi/m 2 . 

Now suppose that at an early stage of the Universe evolution, when the temperature is 
high, say T > m, a gravitating boson ensemble exists in a purely gaseous phase with the 
particle- number density of the order n t ^ ~ m 3 within a volume of large radius -R gas - During 
the cooling down to temperatures T <C m the condensation will take place, in which a 
number of boson stars will be formed. During this evolution the particle number N ~ R 3 as m 3 
is conserved. Since a typical particle number in a boson star is of the order 1/r], the number 
of boson stars formed in the condensation process is approximately n ~ i]N ~ rjR^m 3 
Therefore, the volume occupied by each of the boson stars is of the order R^/n and the 
corresponding radius of the order 

*~%~Htf5-- ( 71 ) 

This quantity is much larger than the radius Rq of the boson star itself, so the condensate 
will occupy only a small portion of the volume V ~ R 3 . At T ~ m, the thermal contribution 
to the particle number is iVth ~ R 3n th ~ I/ 7 ?- Therefore, the contributions to the total 
charge of the thermal part and of the condensate are of the same order of magnitude despite 
the apparent incommensurability of the two densities in 



5 Results and discussion 

To integrate equations (^5|)-(^8l) numerically, we take the radius of the system as R = 
^-i/3 m -i^ j n accorc i ance with (|7ID , and, for definiteness, we take the mass of the boson as 
m = lGeV. The choice of m will practically not affect our results as long asm< Mpi- 

For given a, the initial values Xo an d <fo that yield the ground-state (zero- node) solution 
with the required particle number are not always uniquely determined. The reason is that, 
owing to general-relativistic effects, there exists a finite xo for which iV and M reach a 
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maximum. This limiting configuration is similar to the Oppenheimer-Volkoff limit for the 
degenerate fermion stars. Further increase of Xo makes N oscillate about a nonzero limiting 
value iVoo that corresponds to the infinite central density. In Fig. ^ we plot the particle 
number N versus xo for various a, including a = oo which corresponds to zero temperature. 
It is clear that for each a there will be at least two configurations with the same N, if N is 
slightly below the maximum. However, the configurations represented by the points on the 
right of the maximum on each curve are thermodynamically unstable since their free energy, 
as we shall shortly see, is larger than the free energy of the corresponding configurations on 
the left of the maximum. Therefore, if we restrict our consideration to stable configurations, 
the quantities xo and ipo will be uniquely determined. 

Fixing N = 0.5, which is below the maximum of the zero temperature curve, we now 
plot the initial values Xo an d y?o as functions of inverse a in Figs. [3] and f|, respectively. At 
the critical point l/a c = T c /m = 0.582, xo changes the sign and (p vanishes. Thus, the 
configurations for which a < a c are purely thermal and those with a > a c are mixed. For 
each a we can find u, /i, and T by making use of equations (|54D and 

In Fig. [| we plot the free energy per particle, calculated using (^) and (^), as a function 
of temperature. The full line represents the free energy of the configurations corresponding 
to the points on the left of the maxima in Fig. whereas the dashed line represents the free 
energy corresponding to the points on the right of the maxima. As expected, the latter free 
energy is everywhere larger than the former, hence, the configurations it represents are ther- 
modynamically unstable. In what follows we restrict our attention to stable configurations 
only. The critical temperature T c = 0.582 m is indicated by a vertical line. In Figs. |] and |7| 
the entropy per particle and the total mass are plotted as functions of temperature. Since 
S and M are both continuous functions at T c , we conclude that the process of condensation 
is a second-order phase transition with the naturally defined order parameter ip plotted in 
Fig. |. 

In order to facilitate the physical understanding of our results, let us mention the prop- 
erties of the particular soliton configuration in physical units. The radius of the soliton at 
T = is a few times larger than (cm/h)' 1 , so for m = 1 GeV it is approximately 1.2 x 10~ 12 
cm, whereas the radius of the whole star at T ^ is 1.04 cm. The total mass of the soliton at 
T = is M = 1.3 x 10 11 kg, the particle number N = 7.5 x 10 37 , and the critical temperature 
T c = 0.582 GeV. 

At this stage it is instructive to compare our results with the approximate analytical 
result of a corresponding free-boson system. Taking the lowest contributions in the high- 
temperature expansion, one can find simple analytical expressions that describe properties 
of the free boson gas 0. In particular, the critical temperature is given by 

MIT- 



Using V = AnR 3 /3 we find 
with R measured in units of r]~ 1 ^ 3 m~ 1 and N in units of Within the same approximation 



'< •"> i? -3/2 iV l/2 ) (73) 
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the temperature dependence of the chemical potential is given by 



m \T 



(74) 



for T > T c ^> m, and [l = m for temperatures T < T c . 

In Figs. |8| and |9| we plot the chemical potential /i at infinity and the chemical potential 
p,(0) at the origin, respectively, as functions of temperature. For comparison, in Fig. || we 
also give the discussed approximate analytical result for a free boson system. It is worth 
noting the following two effects. First, the critical temperature is almost exactly equal 
to that obtained for a free relativistic Bose system. Second, the functions describing the 
temperature dependence of \i and /i(0) for T > T c are almost exactly the same. The reason 
for both effects is that goo( r ) in a purely gaseous phase is practically constant on the interval 
(0,i?). The departure of the analytical result for a free boson system from the numerical 
result for a gravitating boson system is a general- relativistic effect below T c , and above T c 
is just due to the high-temperature approximation. 

Our results depicted in Figs. ||-|9| depend, of course, on the choice of R and N. However, 
as in the case of a free gas, the critical temperature should not depend on R nor on N if the 
average density is kept fixed. We have checked that, if we keep the ratio N/R? fixed, our T c 
remains constant to very high accuracy with R varying by several orders of magnitude. 

So far, we have studied the condensation process assuming m <C M P1 , i.e., 77 <d 1 . Now, 
let us investigate an extreme general-relativistic case when the value of the boson mass is 
a sizable fraction of the Planck mass, say m = O.lMpi, yielding rj = 0.01. In this case, the 
natural length scale r q~ 1 / 2 m~ 1 of the purely thermal phase becomes comparable with the 
natural size r/ _1 / 3 m _1 of the mixed phase, and we expect large gravitational effects. 

The calculations are presented in Figs. ID and |TT] in which, respectively, the free energy 
and entropy are shown as functions of temperature for i] = 0.01 and R = lO^^-^m^ 1 . 



The part of the curve in Fig. |10| that starts from T = 0, makes the loop, and ends at the 
sharp cusp, represents the mixed soliton-gas phase. The rest of the curve, which can be 
continuously extended to high temperatures, represents the purely gaseous phase. In the 
temperature interval T = (0.003 — 0.015)m there are three distinct solutions of which only 
two are physical, namely, those for which the free energy assumes a minimum. The cusp 
corresponds to the naive transition point. However, the actual transition takes place at the 
temperature T c , where the free energy of the mixed phase and that of the gas become equal. 
Hence, it is a first-order phase transition. The dotted curves in Figs. |H] and |ll] represent the 
physically unstable solution. In our example, the transition temperature is T c = 0.00607m, 
as indicated in the plots by the dashed line. The latent heat per particle released during the 
phase transition is given by the entropy difference at the point of discontinuity 

AM AS , N 

— — = — T c = 0.0206m. 75 

This phase transition is similar to the gravitational phase transition in self-gravitating 



fermionic systems [21], with one important distinction: in a gravitating bosonic system 



the phase transition is first order only in the extreme general-relativistic regime, whereas in 



gravitating fermionic systems it remains first order even in the Newtonian regime [^, ^3 
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6 Conclusions 



In this work, we have studied a canonical system of self-gravitating bosons at finite tem- 
perature in a general-relativistic framework. We have numerically solved the system of 
self-consistency equations consisting of a Klein-Gordon equation coupled to Einstein field 
equations and the equation of state for a gravitating boson gas. We have investigated the 
circumstances under which this system undergoes a Bose-Einstein condensation. This con- 
densation is quite distinct from the usual one in that it involves the formation of a soliton 
with the spherically symmetric matter distribution concentrated around the origin, as op- 
posed to the usual spatially homogeneous Bose-Einstein condensate. In the T — > limit, 
the soliton becomes a mini boson star. The condensation begins with a second order phase 
transition at the critical temperature T c . The system exists in two phases: a gaseous phase 
above T c and the mixed soliton-gas phase below T c . The critical temperature is approxi- 
mately proportional to the square root of the average particle-number density and is very 
close the the corresponding critical temperature of a free boson gas. 

General relativistic effects become important when the boson mass is a few orders of mag- 
nitude away from the Planck mass. In that case, the condensation begins with a first-order 
phase transition that qualitatively resembles the gravitational phase transition of fermionic 
matter. 

Acknowledgments 

This work was supported by the Ministry of Science and Technology of the Republic of 
Croatia under Contract No. 00980102. 

References 

[1] J.I. Kapusta, Phys. Rev. D 24 (1981) 426. 

[2] H.E. Haber and H.A. Weldon, Phys. Rev. Lett. 46 (1981) 1497. 

[3] H.E. Haber and H. A. Weldon, Phys. Rev. D 25 (1982) 502. 

[4] R. Friedberg, T.D. Lee, and A. Sirlin, Phys. Rev. D 13 (1976) 2739. 

[5] S. Coleman, Nucl. Phys. B 262 (1985) 263. 

[6] D.J. Kaup, Phys. Rev. 172 (1968) 1331. 

[7] R. Friedberg, T.D. Lee and Y. Pang, Phys. Rev. D 35 (1987) 3640. 

[8] P. Jetzer, Phys. Rep. 220 (1992) 163. 

[9] E. Seidel and W.-M. Suen, Phys. Rev. Lett. 72 (1994) 2516. 

[10] E.W. Mielke and F.E. Schunck, Nucl. Phys. B 564 (2000) 185. 



14 



[11] G. Ingrosso and R. Ruffini, Nuovo Cim. B 101 (1988) 369. 
[12] M. Laine and M. Shapshnikov, Nucl. Phys. B 532 (1998) 376. 
[13] G.W. Gibbons and S.W. Hawking Phys. Rev. D 15 (1977) 2752. 

[14] S.W. Hawking, The path-integral approach to quantum gravity, in General relativity: 
An Einstein centenary survey, edited by S.W. Hawking and W. Israel (Cambridge Uni- 
versity Press, Cambridge, 1979). 

[15] J. I. Kapusta, Finite-temperature field theory (Cambridge University Press, Cambridge 
1989). 

[16] R.C. Tolman, Relativity, thermodynamics and cosmology (Clarendon, Oxford, 1934). 

[17] L.D. Landau and E.M. Lifshitz, Fluid mechanics (Pergamon, Oxford, 1959). 

[18] J. Ehlers, Survey of general relativity theory, in Relativity, astrophysics and cosmology, 
edited by W. Israel (D. Reidel, Dordrecht/Boston, 1973). 

[19] J.M. Bardeen, B. Carter and S.W. Hawking, Commun. Math. Phys. 31 (1973) 161. 

[20] N. Bilic and R.D. Viollier, Gen. Rel. Grav. 31 (1999) 1105. 

[21] N. Bilic and R.D. Viollier, Eur. Phys. J. C 11 (1999) 173. 

[22] J. Messer, J. Math. Phys. 22 (1981) 2910. 

[23] N. Bilic and R.D. Viollier, Phys. Lett. B 408 (1997) 75; 

N. Bilic and R.D. Viollier, Nucl. Phys. B (Proc. Suppl.) 66 (1998) 256. 



15 



.n/n(0) 


I 1 11 
- 






\ \ 




\ \ 




\ \ \ 




\ \ 




\ \ 




\ \ 




\ ^ 
\ \ 
\ \ 




\ \ 
\ \ 




\ ^ 




\ \ 




V,„, v , ,; 


In ll 1 . 1 1 


"' 10° 10' 


10 2 



r 



Figure 1: Normalized particle- number density in a boson star at T = for ipo equal to 
0.0764 (full line), 0.01 (long dashed line), 0.001 (short dashed line), and 0.0001 (dotted line). 
The radius r is in units of mr x 




Figure 2: Particle number iV (in units of Mpjm 2 ), versus Xo fc> r various a. 
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Figure 3: Initial value of the metric field xo that yields the zero-node soliton solution with 
iV = 0.5MpJm 2 , within a sphere of radius R = (Mpi/m) 2 / 3 m _1 , versus 1/a. 
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Figure 4: Initial value <fo as a function of 1/a for iV and i? as in Fig. |3|. 
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Figure 5: Free energy per particle F/N (in units of m) versus temperature T (in units of 
for N and R as in Fig. |3]. 




Figure 6: Entropy per particle versus temperature (in units of m) for N and R as in Fig 
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Figure 7: Total mass (in units of M|j/m) versus temperature (in units of m) for N and R 
as in Fig. [|. 




Figure 8: Chemical potential fi (in units of m) at r = oo as a function of temperature (in 
units of m) for gravitating bosons (solid line) and for free bosons in the high-temperature 
approximation (dashed line) for N and R as in Fig. || 
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Figure 9: Chemical potential /i (in units of m) at r = as a function of temperature (in 
units of m) for N and R as in Fig. |3|. 




Figure 10: Free energy per particle (in units of m) versus temperature (in units of m) for 
r] = m 2 /M^ = 0.01 and N = 0.5M^/m 2 . 
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Figure 11: Entropy per particle versus temperature (in units of m) for r\ = m 2 /Mpj = 0.01 
and N = 0.5M^/m 2 . 
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